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Interaction between shrub encroachment and water 


infiltration on a hillslope at the typical Steppe 


Abstract The interaction between the surface hydrologic cycle and the 
shrub-encroached landscape at different slope positions remains poorly 
investigated. This study aims to explore the interaction between the water 
infiltration patterns affected by shrub encroachment at different hillslope 
positions. Soil water content and temperature were continuously measured 
at 10-min intervals at four or five depths under shrub patches and the grass 
matrix at four slope positions of a Caragana microphylla encroached 
hillslope from July 2009 to May 2013. The rainfall and meltwater 
infiltrations were estimated based on above data. Results showed that the 
rainfall infiltration ratios (IRs) at the grass matrix were as high as 
0.78+0.08, except at the lower site, where it was only 0.47. The IRs of 
shrub patches increased from 0.38 at the top site to 0.77 at the lower site. 
The IRs were higher at the grass matrix than that at the shrub downslope 
edges at the top, upper, and middle sites of the hillslope due to the raised 


microtopography of the shrub mounds. However, at the lower site, IR was 
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higher at the shrub patch than that at the grass matrix than due to more 
upper slope runoff input and higher infiltration capacity at the shrub 
patches. The preferential flow was not an important factor influencing the 
redistribution of water resources on the slope. Snow and ice were blown 
up by wind and accumulated in the shrub patches and their lees resulted 
higher water input to the shrub patches than that in the grass matrix during 
snowy years. Shrub encroachment changed the microtopography, soil 
property under different canopy and slope positions, and further affected 
the surface hydrological processes. The feedbacks between shrub 
encroachment and water infiltration varied at different sites of the hillslope 


and affected the development of shrub patches. 


Key words: Infiltration; rainfall; meltwater; shrub encroachment 


1 Introduction 

Shrub encroachment has been a worldwide phenomenon over more 
than 150 years in the arid or semi-arid regions, which refers to the increase 
in plant density, coverage and biomass of shrub in the grassland ecosystem 
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(Van Auken, 2000; Gibbens ef al., 2005; Eldridge et al., 2011; D'Odorico 
et al., 2012). It is believed that erosion and climate change are the main 
causes of shrub encroachment (Sturm et al., 2001a; D'Odorico et al., 2012). 
Shrub encroachment redistributes light, heat, water, nutrients, and other 
resources, and changes net primary production, species diversity, etc. 
(Eldridge et al., 2011; Peng et al., 2014; Li et al., 2019; Zhou et al., 2019). 
The encroached vegetation interacts strongly with the environment (Rolo 
and Moreno, 2019; Turnbull and Wainwright, 2019), and their feedbacks 
would promote the encroachment and result a new longtime stable “climax 
community” (Thomas and Pittillo, 1987). The characteristics of the shrub 
patches were related to the soil water content and local precipitation, and 
with the increase of water resources, the optimal shrub coverage also 
increased (Valentin and D'Herbés, 1999; Zheng et al., 2015a). 

Positive interactions between shrub encroachment and water resource 
of shrub patches were widely found. Li et al. (2013a) reported that soil 
water content at the depth of 60 cm under shrub patches was higher than 
that under grass matrix, which was related to their different CaCO; 
horizons characteristics. Higher infiltration depths and more preferential 


5 / 48 


75 


76 


77 


78 


79 


80 


81 


82 


83 


84 


85 


86 


87 


88 


89 


90 


91 


92 


flows under shrub patches were found in dying experiments with the same 
volume of dye solution (Zhang et al., 2012; Peng et al., 2013), which also 
revealed higher infiltration capacity under shrub patches (Peng et al., 2013; 
Li et al., 20132). Shrub patches usually had lower runoff rates than grass 
land (Pierson et al., 2010; Munoz-Robles et al., 2011; Peng et al., 2013). 
Soliveres and Eldridge (2014) and Marquart et al. (2019) also reported that 
steady-state infiltration and sorptivity were greater under shrubs canopies 
than those in open areas, and positive effects of shrub were ascribed to the 
reduction of trampling from livestock. Hu et al. (2015, 2018) found that 
soil macroporosity and macropore volume in shrubs patches were higher 
than those in grass patches, and they increased with increasing shrub root 
network density. Thus, it was speculated that water could transport to the 
deeper soil layers in shrub patches (Hu et al., 2015; 2018). The root 
channels of the shrub were thought to be preferential pathways of water 
movement and could let water reach deep soil layer (Martinez-Meza and 
Whitford, 1996; Wang et al., 2007; Li et al., 2009; Zhang et al., 2012). On 
the other side, Rossi and Ares (2017) modeled the water fluxes between 
inter-patches and shrub mounds and found that shrub mounds could 
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directly obtain more water resource from the rain, but overland flow from 
the inter-patches could only reach the border of the shrub mounds and 
could not route onto the shrub mounds, and no ponded-water inundates the 
mound. So, as most of the former infiltration and runoff studies were 
conducted under artificial condition, how the infiltration in shrub patches 
act under natural precipitation, and whether the effects are positive in 
different parts of the shrub patches (center and edge) as well as on different 
positions of a hillslope were still unclear. 

The snow can be an important soil water input and influence soil water 
infiltration and storage (Stewart et al., 2004; Murray and Buttle, 2005). 
Meltwater can form macropore flow in rooted soil, which effectively 
storages water for plants (Newman et al., 2004). Snow melting in arid and 
semi-arid areas can improve soil moisture, create conditions for seed 
germination and growth, and even have a profound impact on the water 
competition within and between species and on the spatial pattern of 
vegetation (Li et al., 2006). The redistribution of water resource will affect 
the productivity and distribution of vegetation (Weltzin et al., 2003). In the 
Arctic, deciduous shrubs obtained more water from snow than grass in the 
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growing season, and this meltwater-focused feedback loop could explain 
the shrub expansion (Jespersen et al. 2018). In addition, shrub 
encroachments have greatly changed the landscape pattern of grassland 
ecosystem, local radiation balance, physical characteristics and structures 
of snow cover (Murray and Buttle, 2003), and further changed the local 
microclimate, runoff, and soil water conditions (Wilcox et al., 2003; Li et 
al., 2006). The biomass and height of shrub patches were higher than those 
of the grass inter-patches in seriously disturbed conditions (Peng et al., 
2013). Thus, shrub patches could trap and retain more snow and meltwater 
resource (Sturm et al., 2001b; Tape et al., 2006; Myers-Smith et al., 2011; 
Yan et al., 2019). But how these effects perform on different slope position 
and under different snow level in the north China need more investigation. 

Therefore, the aims of this study were to (a) investigate the infiltration 
characteristics under natural rainfall at different sites of a slope with shrub 
encroachment, and (b) evaluate the effects of the meltwater infiltration on 
water resources of shrub patches. This study would explore the 
relationships between the development of the encroached shrub patches 
and the infiltration pattern at different sites of the slope. The results are of 
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great significance to explain the spatial distribution characteristics of 
encroached shrub and predict its future development trend under changing 


climate background. 


2 Methods and Materials 
2.1 Study site 

This study was carried out at the Taipus Banner, Inner Mongolia, China. 
This area belongs to the typical Eurasian Steppe and has continental 
temperate semiarid climate. Its mean annual temperature and mean annual 
precipitation are 1.6 ?C and 392 mm, respectively (Li et al., 2013a). The 
hottest month is July (17.8 ?C on average) and the coldest month is January 
(-17.6 °C on average). 65% of precipitation occurs from July to September. 
Annual pan evaporation reaches 1900 mm on average. The average annual 
wind velocity is 3 - 5 ms! and the annual period of force-seven wind and 
above is 20 - 80 d, mainly occurring in spring and winter. The annual frost- 
free period is around 100 d, from late May or early June to early September. 
The zonal soils are chestnut soil and light chestnut soil, which are 
equivalent to Calcicorthic Aridisol according to USDA Soil Taxonomy. 
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Soil textures were mainly sand and sandy loam. The dominant plant species 
were Stipa krylovii Roshev., Cleistogenes squarrosa (Trin. Keng, 
Artemisia frigida Willd., and Leymus chinensis (Trin.) Tzvel (Qiu et al., 
2011). Caragana microphylla Lam expansion increased during recent 
decades in this region. It is estimated the area of C. microphylla encroached 
grassland amounted to more than 5.1 x 10° ha (Zhang et al., 2006), and the 
coverage of this shrub can reach 40.1% (Peng et al., 2013). 

One C. microphylla encroached hillslope was selected for this study 
(Figure 1). This west slope is 128 m in length and has an average slope 
gradient of 12.5% (Li et al., 2013a). This slope was evenly divided into 
four sites from the ridge to the valley: the top, upper, middle, and lower 
sites (Figure 1). The sizes of shrub patches were 3.93, 5.86, 5.84, and 7.68 
m? at the top, upper, middle, and lower sites of the slope, respectively, and 
their canopy heights were 0.43, 0.55, 0.47, and 0.41 m, respectively (Li et 
al., 20132). The above ground biomasses of the shrub patches at four sites 
were 343, 400, 403, and 428 g m? from top to lower sites, respectively, 
while that of the grass matrix were 108, 140, 141, and 185 g m?, 
respectively (Li et al., 20132). The canopy coverage of shrub patches was 
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93-95%, and that of the grass matrix is 68-75% in the peak season (Li et 
al., 2013a). The rooting depths were 60 - 70 cm and 20 - 25 cm at the shrub 
patches and at the grass matrix, respectively (Li et al., 2013a). 

«Figure 1> 
2.2 Soil water content and temperature measurement 

Eight soil profiles were instrumented at the top, upper, middle, and 
lower sites of the slope, two for each site (Li et al., 2013a): one at the 
downslope edge of a C. microphylla shrub patch, and the other at the grass 
matrix, about 1 m from the upslope edge of the C. microphylla shrub patch 
(Figure 1). Five ECH2O STE sensors (Decagon Devices, Pullman, 
Washington, USA) were installed in each soil profile to monitor soil water 
content (0, m? m^?) and soil temperature (T;, °C) at the soil depths of 10, 
20, 40, 60, and 100 cm. Except for the middle site, the ECH2O STE sensors 
was replaced by Decagon Drain Gauges at the100 cm depth to detect soil 
water drainage (Figure 1). A total of thirty-eight ECH2O STE sensors were 
installed at four sites of the slope. All data of 0 and T;, as well as drainage, 
if available, at each profile were logged into a 5-channel Decagon's Em50 
data logger at 10-min intervals. These ECH2O STE sensors can provide 
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accurate, precise and continuous measurements, with the accuracy of 
+0.012 m? m? for moisture, +0.3°C for soil temperature (Czarnomski et al., 
2005; Li et al., 20132). The experiments started on June 8, 2009, and data 
since July 1, 2009, or later was used considering the disturbance of 
profiling and setting of the sensors (Table 1). The monitoring last to May 
16, 2011, at the top, upper, and lower sites, while at the middle site, it last 


to May 21, 2013 (Table 1). 


«Table 1> 


2.3 Meteorological observations 

A Dynamet weather station (Dynamax Inc., Houston, TX, USA) was 
set up to measure the meteorological data at the center of the slope. Rainfall 
was measured using a tipping bucket rain gauge (Model TE525, Campbell 
Scientific, Logan, UT, USA). Air temperature (7;, °C) was recorded using 
a temperature and humidity probe (Model CS500, Campbell Scientific, 
Logan, UT, USA) at 2 m above the ground. Wind speed and direction were 


detected by a windvane and anemovane (Model 03001, RM Young Co., 
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Traverse City, MI, USA), positioned at 3 m above the ground. Snowfall 
data were obtained from the Weather Station of Taipus Banner (41.88°N, 
115.27°E,1468.9 m a.s.l.). 
2.4 Infiltration ratios estimation 

The rain events were discretized by assuming a minimum interevent 
time of 24 h. Soil water storage (SWS, mm) of a 1 m profile was calculated 
as: 

SWS = $0,D, (1) 

where 0 (m? m?) represents 0 measured by 5TE seniors at soil layer k, D; 
is the depth of 0 represents (mm). n represents the number of the STE 
sensors at the 1 m profile. At the top, upper, and lower sites, n = 5, and D; 
= 150, 150, 200, 300, and 200 mm for k = 1, 2, 3, 4, and 5, respectively; 
while at the middle sites, n = 4, and D, = 150, 150, 200, and 500 mm for k 
— 1, 2, 3, and 4, respectively. Infiltration amounts (/, mm) in a rainfall event 
was estimated as the maximum difference of SWS before and after the 
rainfall event (before the next rainfall). The infiltration ratio (ZR) was 
estimated as the slopes of the linear regression of / against rainfall amounts 
at the experimental period (Zhang and Li, 2018). 
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2.5 Preferential flow detection 

The preferential flow (PF) during natural precipitation were detected 
by the response sequences of 0 of different soil layers. Only events with 
over | mm of rainfall were analyzed considering that small rainfall could 
not infiltrate deeper than the first sensor. The 0 was examined and thought 
tohave response to rainfall if its increase was over 0.01 m? m? during an 
rainfall event due to the precision and noise levels of the sensors was 0.01 
m? m? (Graham and Lin, 2011). If a deeper sensor responds before, or 
simultaneously with a shallower sensor, or when a shallower sensor does 
not respond, it was thought that preferential flow happened at this profile 
in this event (Graham and Lin, 2011; Hopkins et al., 2016). 
2.6 Meltwater infiltration estimation 

When the soil was frozen, 0 decreased as soil water transferred into ice. 
During the frozen period, soil water was protected from evaporation. After 
thawed, 0 can nearly return to the same level before soil was frozen. Basing 
on this assumption and ignoring the evaporation during this period, 
meltwater infiltration was calculated as the differences of the maximum 
SWS before the first rain in a spring (after thawing) and the SWS just before 
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the soil freezing in the last winter (Figure 2). 


«Figure 2> 


3 Results 
3.1 Precipitation and wind characteristics 

Characteristics of rainfalls in the experimental area during the nearly 
four years (from July 1, 2009, to May 21, 2013) were given in Figure 3. A 
total of 144 rainfall events (excluding snowfall) occurred in nearly four 
years, accumulating 904.71 mm of rainfall. They were mainly short small 
storms, with 46 (31.94%) rainfall events smaller than 1 mm, and 96 
(66.67%) rainfall events not higher than 5 mm. Rainfalls larger than 10 mm 
only accounted for 13.19% of all rainfall events but accounted for 65.52% 
of the rainfall amount. The maximum rainfall was 59.69 mm recorded 
during September 15 — 18, 2010. The rainfall intensity ranged from 0.02 to 
13.20 mm h'!. The maximum 10 - min intensity ranged from 0.20 to 73.15 
mm h, with only 10.42% of which larger than 10 mm h'!, and only one of 
them was larger than the 63.23 mm h“ infiltration rate of grass matrix at 
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suctions of 0.5 mm (Li et al., 2013a). The snowfalls were 39.5, 49.5, 29.1, 
and 45.1 mm during the winter and spring of 2009-2010, 2010-2011, 2011- 
2012, and 2012-2013, respectively. The wind in winter and spring were 
mainly from W, WSW, SW, and SSW, and their speeds mainly ranged from 
2 - 10 m s with an average speed of 5.2 m stand a maximum of 18.8 m s 


! (Figure 4). 


«Figure 3> 


«Figure 4> 


3.2 Rainfall infiltrations 

Rainfall infiltrations varied at different slope position and under 
different land covers. At the top, upper, and middle grass matrix, /Rs were 
as high as 0.7840.08, while at the lower grass matrix, it was only 0.47 
(Figure 5). The /Rs of shrub patches increased from 0.38 at the top site to 
0.77 at the lower site (Figure 5). 

The yearly total infiltration showed similar results to the /Rs (Figure 6). 
At the top and upper sites, the grass matrix had higher yearly total 
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infiltration than shrub patches, with higher differences at the top site than 
at the upper site. The yearly total infiltration at the lower shrub patch was 
higher than that at the lower grass matrix. 

The yearly infiltrations in 2010 were available at all sites with a 
minimum of 69.4 mm at the top shrub patch and a maximum of 194.5 mm 
at the upper grass matrix. In 2011 and 2012, the whole year infiltrations 
were available at the middle sites. The yearly infiltration of2011 was 212.4 
mm at the middle shrub patch, the maximum of all sites in the whole 
experimental period, and much higher than the 138.3 mm at the middle 
grass matrix. In 2012, the yearly infiltration at the middle shrub patch was 
123.1 mm, lower than that of 155.6 mm at the middle grass matrix. From 
the 2010 to 2012, the yearly infiltration at the middle shrub patch was 152.1 
+52.2 mm, insignificantly higher than that at the middle grass matrix, 


143.4+ 10.5 mm, basing on the 2-tailed paired sample T test (p=0.818). 


«Figure 5> 


«Figure 6> 
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PF ratios were only 4.1% and 2.4% at grass matrix and shrub patches, 
respectively (Table 1) with total 19 PF events recorded. Most of PF events 
were recorded at 60 and 40 cm at the grass matrix and shrub patches, 
respectively (Table 1). Grass matrix had higher PF ratios than the shrub 
patches. There were no PF events recorded at the top sites. At the grass 
matrix, the most PF events were recorded at the middle site, followed by 
the lower site; while at the shrub patches, the highest PF ratios were 
recorded at the middle site. All these 19 PF events came from different 
rainfall events, with the ranges of rainfall depth, rain intensity, and the 
maximum 10 - min rainfall intensities of 1 - 37.85 mm, 0.05 - 9.60 mm hr 
1 and 1.20 - 57.90 mm h', respectively. 

3.3 Meltwater infiltrations 

The meltwater infiltrations varied in different years (Figure 7). In the 
winters of snowless years (2009 - 2010 and 2011 - 2012), the meltwater 
infiltration was low, no more than 5 mm except at the middle grass matrix 
site (6.3 mm), and no infiltration at the shrub patches at both the upper and 
lower sites. In these years, the meltwater infiltration was higher at the grass 
matrix than that at the shrub patches. While at the snowy years, the 
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meltwater infiltration was relatively higher, ranging from 5.17 to 12.94 mm 
in the spring of 2011 (Figure 7). In these snowy years, the meltwater 
infiltration was higher in the shrub patches than that in the grass matrix, 
except for the middle site in the winter of 2010-2011. 

Especially, meltwater infiltration was quite high at the middle grass 
matrix and shrub patch in the spring of 2013 (Figure 8). The soil 
temperature dropped chronically below 0 °C from 1 a.m. and 6 a.m., 
November 13, 2012 at the middle grass matrix and the shrub patch, 
respectively. Before the soil froze in the middle grass and shrub sites, 
their SWSs were 137.6 and 102.91, respectively. After the soil was frozen, 
no rainfalls were recorded. As soil water transferred to ice, the liquid 0 
successively reduced from 10 cm to 60 cm depth , and the liquid SWS of 
the whole profile gradually reduced. The positive soil temperature at 10 cm 
depth was first recorded at 4 p.m. March 7, 2013, at the middle grass matrix, 
while at the middle shrub patch, it was at 12 a.m. March 8, 2013. During 
these two days, the air temperature increased rapidly and reached 12.0 °C 
at 3 p.m. After the soil temperature at 10 cm depth stride over the freezing 
point, the 0 increased quickly and successively at each soil layer. The liquid 
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SWS of the whole profile bulged rapidly at both middle grass and shrub 
sites, reaching 224.4 and 232.7 mm at 7 p.m. March 8, 2013, respectively. 
The meltwater infiltration reached 86.8 mm and 129.8 mm at the middle 


grass and shrub sites, respectively. 


<Figure 7> 


<Figure 8> 


4 Discussions 

Lower /Rs and PF ratios were found at the downslope edge of shrub 
patches than that at the grass matrix except for the lower site. This can be 
related to the slope topography and position where the STE sensors were 
installed. The shrub patches form mounds and are raised microtopography 
(Figure 1 and Figure 9). This is because that shrub prevents the ground 
from livestock trampling and erosion, and on other side accumulates wind 
and water sediment (Stavi et al., 2008; Lu et al., 2019; Yang et al., 2019). 
In this condition, the water flow from upslope would first reach the upslope 
edge of shrub mounds, and water flow would mainly infiltrate there or 
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bypass the of shrub mounds, which results less water source for the 
downslope edge of the shrub patches. The 5TE sensors were installed at 
the downslope edge of the shrub patches and they had the relatively higher 
slope gradients (Stavi et al., 2008), which had higher possibility to generate 
more runoff and reduce water infiltration. Therefore, the microtopgraphy 
of shrub patches resulted that the downslope edge of the shrub patches 
generated more local runoff and stored less water source from upslope. 
This would be the main reason for the low ZR there, which is consistent 
with the result of Li et al. (2013b), who reported lower soil moisture at the 
downslope edges of shrub patches comparing to that at the upslope edges 
after rainfalls. 

Higher JR was found at the downslope edge of shrub patches at the 
lower site. This is because there was more runoff from the upslopes and 
there was plenty of water flowing over or around the shrub mound to gather 
at the downslope edge of the shrub patches. The upslope runoff would form 
ponding water and increase their infiltration as the shrub patches had higher 
infiltration rates (Li et al., 2013a) and deeper infiltration depth under water 
ponding conditions (Zhang et al., 2012; Peng et al., 2013). In addition, the 
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shrub patch at the lower site had the highest aboveground biomass and root 
depth, higher soil organic matter (Li et al., 2013a), and more sand and silt 
(Zhang et al., 2019), which could generate more macropores and increase 
water infiltration (Hu et al., 2015; 2018). However, the lower grass matrix 
site had the lowest total infiltration among the four grass matrix sites. This 
can be explained by more clay deposited on the lower grass matrix, which 
could block the macropores and reduce infiltration capacity (Li et al., 
2013a). 

The PF ratios at the experimental sites were no more than 5%, much 
lower than the 17 to 54% at the Shale Hills Critical Zone reported by 
Graham and Lin (2011). This may be related to their different soil types 
and rainfall characteristics. The soil of Shale Hills Critical Zone was 
formed by shale residuum and colluvium, which were not very well 
developed due to the cold weather and contained a lot of rock which 
increased the formation of PF. While the soil at our experimental sites was 
chestnut soil and light chestnut soil, developed from aeolian sediment, and 
had relative uniform texture throughout the profile. This would not be 
beneficial to the formation of PF. Another reason is that there were few 
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heavy rainfall events in the experimental sites. The most influential factors 
of PF were water input amounts and peak input intensity (Hopkins et al., 
2016). The PF usually occurrences when the water input over a threshold 
(Mcgrath et al., 2010). As showed in Figure 3, rainfalls larger than 10 mm 
only accounted for 13.19% of all rainfall events. The more aboveground 
biomass at the shrub patches(Peng et al., 2013) also intercepted more 
rainfall and reduced the net rainfall reached the ground (Zhang et al., 2018). 
The lower PF ratios of shrub patches can be connected to their low /Rs, 
which equaled to their lower water input. 

Meltwater infiltration is also from two sources: One is the snow and ice 
on the site, and the other is the runoff from the upslope. In the snowless 
years, the meltwater flow was weak and bypassed the shrub mounds. Thus, 
the grass matrix could get meltwater from the upslope while the downslope 
edge of shrub mounds could not get meltwater, resulting in lower meltwater 
infiltration at the downslope edge of shrub compared with the grass matrix. 
So, the meltwater infiltration was lower at the downslope edge of shrub 
than at the grass matrix. On the contract, in the snowy years, snow and ices 
were trapped and retained on the shrub patches and at the lee side. The 
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strong meltwater flow from the upslope and the top of the shrub mound 
could reach the downslope edge of the shrub patch and then increase the 
meltwater infiltration there (Figure 9). Due to the higher infiltration 
capacity of the shrub patches, they had higher meltwater infiltration. In 
general years, the meltwater infiltration did not exceed 15 mm, which could 
only be about 20 mm of rainfall infiltration. However, in some separate 
years, snow was heavy and the shrub patches retained a lot of snow and ice. 
Additionally, if the weather in spring warm up rapidly, the snow and ice 
stack would melt rapidly and generate flood over the slope (Figure 9). 
Compared with most years, this meltwater flow could infiltrate more 
(Figure 7) and reach the soil layer of 60 cm or deeper (Figure 8). 

The meltwater infiltration increased the 0 in the spring. In the snowless 
years, Jonly weakly increased and decreased quickly below the level 
before the soil was frozen in the last year. The effect of the meltwater, the 
period when 0 in spring (before the first rain) is higher than 0 before frozen 
in the last winter, was usually less than 10 days. However, at the snowy 
years, Ostrongly increased and the effect of 0 was usually longer than 20 
days and lasted to the start of vegetation greenup. Especially, in the spring 
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of 2013, the 0 at the middle grass and shrub sites was higher than that 
before frozen till the end ofthe experiment, lasting more than 73 days. The 
increased 0 of meltwater water infiltration was an important soil water 
resource in the semiarid region, which had a great influence on the sprout 
and regreen of vegetation there (Stewart et al., 2004; Barrere et al., 2018). 
The heavy and deep meltwater infiltration in the spring of 2013 would even 
have a significant impact on the growth of local vegetation (Zheng et al., 


2015b). 


«Figure 9> 


5 Conclusions 

Shrub patches could retain more water resources at the lower site of the 
hillslope and meltwater can be an important water resource of the shrub 
patches during snowy years. The preferential flow might not be an 
important factor influencing the redistribution of water resources on the 
slope. The downslope edge of the shrub patches retained lower rain water 
resources than the grass matrix at the top, upper, and middle sites. The 
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downslope edge of shrub patches at the middle, upper, and top sites would 
degrade due to less water resource, and the upslope edge of shrub patches 
would expand toward the top of the slope. However, as shrub patches 
expand toward the top, the runoff resources and their /Rs declined. Thus, 
the size of the shrub patch would stay steady when the degradation and 
expansion are balanced. The shrub patches at the lower site would expand 
in all direction and expand to larger sizes as they could retain more water 


resources. 
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615 Table: 


616 Table 1 Preferential flow (PF) characteristics of different sites. 


Period Rain events PF events PF ratio (%) 
id Begin End Days Total »1mm 20cm 40cm 60cm 100cm Total 20cm 40cm 60cm 100cm Total 
Top grass Jul. 1,2009 May 16,2011 684 76 54 0 0 0 0 0 0.0 0.0 0.0 0.0 0.0 
Top shrub Aug.1,2009 May 16,2011 6053 70 49 0 0 0 0 0 0.0 0.0 0.0 0.0 0.0 
Upper grass Jul. 1,2009 May 16,2011 684 76 54 1 0 0 0 1 1.9 0.0 0.0 0.0 1.9 
Upper shrub Jul. 17,2009 May 16,2011 668 75 53 0 1 1 0 2 0.0 1.9 1.9 0.0 3.8 
Middle grass Jul. 1,2009 May 21,2013 1420 163 109 1 3 3 - 7 0.9 2.8 2.8 - 6.4 
Middle shrub Jul. 1,2009 = May 21,2013 1420 163 109 1 1 2 - 4 0.9 0.9 1.8 - 3.7 
Lower grass Aug. 1,2009 May 16,2011 653 70 49 0 0 4 0 4 0.0 0.0 8.2 0.0 8.2 
Lower shrub* Jul. 1,2009 May 16,2011 543 70 51 0 1 0 0 1 0.0 2.0 0.0 0.0 2.0 
Grass 2 3 7 0 12 0.7 0.7 2.7 0.0 4.1 
Shrub 1 3 3 0 7 0.2 1.2 0.9 0.0 2.4 


617  F*missing period: Nov.1, 2009 - Mar. 21, 2010 


618 
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Figures: 

Figure 1 Diagrammatic sketch of the experimental hillslope with C. 
microphylla patches interspersed in sparse grass matrix, and sites for soil 
moisture sensors and meteorological tower instrumentation. 

Figure 2 Meltwater infiltration estimation. Soil water storage (SWS, back 
line), and soil temperature (red line) crossing soil frozen and the first rain 
after soil thawed were given. The red dash line shows 0 °C. Black dot lines 
show the frozen soil period. The green dot line shows the SWS just before 
soil frozen and thawed. The blue dot line shows the maximum of SWS after 
soil thawed and before the first rain. 

Figure 3 Characteristics of rainfalls at the experimental area from July 1, 
2009 to May 21, 2013. The grey solid line shows the density of rainfall 
amount and the grey vertical lines show the occurs of rainfall events; 
yellow columns show the count of rainfall events; the blue dash line shows 
the cumulative percent of rainfall amount; the blue bubbles show the 


rainfall amount, rainfall duration, and their sizes show rainfall intensity. 


37 / 48 


635 


636 


637 


638 


639 


640 


641 


642 


643 


644 


645 


646 


647 


648 


649 


650 


651 


652 


Figure 4 Wind speeds and direction in January, February, March, 
November, and December of the experimental period at 3 m above the 
middle of the slope. 

Figure 5 The relationships between the changes of soil water storage 
(ASWS) and rainfall amounts (P). 

Figure 6 Yearly rainfall infiltration at different sites in 2009 and 2010. 
Rainfall infiltration in 2009 were since July 1 except the top shrub site and 
the lower grass site, which were since August 1, 2009. 

Figure 7 The meltwater infiltration at different sites and periods. 

Figure 8 Meltwater infiltration at the middle site in the period of 2012 - 
2013. (a) Rainfall and air temperature (7a) at 2 m above ground; (b) soil 
water storage (SWS) and soil temperature at 10 cm depth at the middle 
shrub site; (c) soil water content (@) at 10, 20, 40, and 60 cm depth at the 
middle shrub site; (d) soil water storage (SWS) and soil temperature at 10 
cm depth at the middle grass site; (c) soil water content (0) at 10, 20, 40, 
and 60 cm depth at the middle grass site. 

Figure 9 Photo of meltwater flow on the slope. Photo taken on March 7, 


2011. Snow and ice were retained at the lee side of the shrub patch and 
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653 melted at the south first. Note that there were snow and ice on the shrub 
654 mound earlier than the time of this photo taken, and they melt earlier due 
655  tothe stems of shrub had lower albedo and could absorbed more radiation. 


656 Photo courtesy of Si-Yi Zhang. 
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Figure 1 Diagrammatic sketch of the experimental hillslope with C. 
microphylla patches interspersed in sparse grass matrix, and sites for soil 


moisture sensors and meteorological tower instrumentation. 
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Figure 2 Meltwater infiltration estimation. Soil water storage (SWS, back 
line), and soil temperature (red line) crossing soil frozen and the first rain 
after soil thawed were given. The red dash line shows 0 °C. Black dot lines 
show the frozen soil period. The green dot line shows the SWS just before 
soil frozen and thawed. The blue dot line shows the maximum of SWS after 


soil thawed and before the first rain. 
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Figure 3 Characteristics of rainfalls at the experimental area from July 1, 
2009 to May 21, 2013. The grey solid line shows the density of rainfall 
amount and the grey vertical lines show the occurs of rainfall events; 
yellow columns show the count of rainfall events; the blue dash line shows 
the cumulative percent of rainfall amount; the blue bubbles show the 


rainfall amount, rainfall duration, and their sizes show rainfall intensity. 
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680 Figure 5 The relationships between the changes of soil water storage 
681  (ASWS) and rainfall amounts (P). 
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683 Figure 6 Yearly rainfall infiltration at different sites in 2009 and 2010. 
684 Rainfall infiltration in 2009 were since July 1 except the top shrub site and 


685 the lower grass site, which were since August 1, 2009. 
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687 Figure 7 The meltwater infiltration at different sites and periods. 
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Figure 8 Meltwater infiltration at the middle site in the period of 2012 - 
2013. (a) Rainfall and air temperature (7;) at 2 m above ground; (b) soil 
water storage (SWS) and soil temperature at 10 cm depth at the middle 
shrub site; (c) soil water content (0) at 10, 20, 40, and 60 cm depth at the 
middle shrub site; (d) soil water storage (SWS) and soil temperature at 10 
cm depth at the middle grass site; (c) soil water content (0) at 10, 20, 40, 


and 60 cm depth at the middle grass site. 
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Figure 9 Photo of meltwater flow on the slope. Photo taken on March 7, 
2011. Snow and ice were retained at the lee side of the shrub patch and 
melted at the south first. Note that there were snow and ice on the shrub 
mound earlier than the time of this photo taken, and they melt earlier due 
to the stems of shrub had lower albedo and could absorbed more radiation. 


Photo courtesy of Si-Yi Zhang. 
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